Anomalous fermion number nonconservation: Paradoxes in the 

level crossing picture 

Y. BurnieiQ 

Institut de Theorie des Phenomenes Physiques, 
Ecole Polytechnique Federale de Lausanne, CH-1015 Lausanne, Switzerland 

(Dated: September 4, 2006) 

Abstract 

In theories with anomalous fermion number nonconservation, the level crossing picture is consid- 
ered a faithful representation of the fermionic quantum number variation. It represents each created 
fermion by an energy level that crosses the zero-energy line from below. If several fermions of var- 
ious masses are created, the level crossing picture contains several levels that cross the zero-energy 
line and cross each other. However, we know from quantum mechanics that the corresponding 
levels cannot cross if the different fermions are mixed via some interaction potential. The simulta- 
neous application of these two requirements on the level behavior leads to paradoxes. For instance, 
a naive interpretation of the resulting level crossing picture gives rise to charge nonconservation. 
In this paper, we resolve this paradox by a precise calculation of the transition probability, and 
discuss what are the implications for the electroweak theory. In particular, the nonperturbative 
transition probability is higher if top quarks are present in the initial state. 
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I. INTRODUCTION 



When a classical conservation law is broken by quantum corrections, It is said that the 
associated symmetry is anomalous. An anomaly in a current associated with gauge sym- 
metry ruins the consistency of the theory. The requirement that all gauge anomalies cancel 
strongly restricts the possible physical theories. On the other hand, anomalies arising in 
other type of currents can lead to interesting physics. For instance in strong interactions, 
the anomaly in the chiral current is important in the well-known pion decay to two pho- 
tons. In weak interactions, there is an anomaly in the baryon number current. Although 
anomalous baryon number violating transitions are strongly suppressed at small energies, 
they could be at the origin of the baryon asymmetry of the universe. 

Anomalous transitions leading to fermion number nonconservation arise in the elec- 
troweak theory or any other model with a similar vacuum structure. The crucial feature is 
the existence of an infinite number of degenerate vacua, separated by energy barriers and 
the transition between them leads to the creation, or destruction, of fermions. The energy 
barrier can be passed by either by tunneling, which is represented by an instanton or 
yy thermal excitations P, 0]. In the second case, the relevant configuration is the sphaleron 
^. It is defined as the maximum energy configuration along the path of minimal energy 
connecting two neighboring vacua. 

To visualize anomalous fermion number nonconservation, let us consider the path in the 
bosonic field space, parametrized by r, that relates two neighboring vacua via the sphaleron 
configuration. If the bosonic fields evolve very slowly along this path, the fermionic states 
can be found by solving the static Dirac equation Hj^^n = En'^n- This equation has positive 
as well as negative energy states. A way to represent the fermionic vacuum state is the Dirac 
sea. All states with negative energy are filled, whereas all positive energy states are empty. 
We are interested in the variation of the Dirac function of r. On a graph containing 

all energy levels as function of r, it may happen that an initially negative (therefore occupied) 
energy level crosses the zero energy line and becomes a real particle. This is the level crossing 
picture representation of the anomalous fermion number nonconservation . In the case of 
the electroweak theory, one level for each existing fermionic doublet crosses the zero-energy 
line in the transition between two adjacent vacua ^. 

The level crossing picture can be thought of as a quantum mechanical description of 
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FIG. 1: Naive picture of level crossing containing only the two levels that crosses the zero-energy 
line. The heavy particle is represented with a thick line. Two cases are pictured; without mixing 
(a) and with mixing between the two fermions (b). In the latter case the energy levels can't cross, 
therefore the heavy particle becomes a light one and vice-versa. 

fermion creation. This description is assumed to match the complete quantum field theory 
when the background evolves very slowly. 

Consider now the case of two fermions where i = 1,2 is the fiavor index. We first 
assume that the different fiavors are not mixed by any interaction term, that is to say the 
Dirac equation, which generally reads Hij'^j = E'^j, can be diagonalized in fiavor space for 
any r. We will call fermions for which H is diagonal, independent. On the level crossing 
picture for two fermions with different masses^, we see that two energy levels cross the zero- 
energy line and cross each other. A simplified level crossing picture containing these two 
levels is given in Fig. ^a-^ 

On the other hand, if the two fermions are mixed, for instance by the interaction between 
them and the background sphaleron fields, and the Dirac Hamiltonian is not diagonalizable, 
we know from quantum mechanics that the energy levels cannot cross each other. Therefore, 
in this case, the heavy fermion becomes the light one and the light one becomes the heavy 
one, see Fig. [T]b. Indeed, generally there are excited states of the light particle and the light 
fermion evolves to one of them, see Fig. |21 Therefore, we conclude from the level crossing 
picture that two light fermions are created in the case with mixing instead of a light and a 
heavy one. 

^ We consider here fermions made massive through their Yukawa coupling to the Higgs field and not by a 
tree mass term. 

^ The full level crossing picture of the theory we consider in the following is given in Fig. [3 
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Suppose now that we introduce another gauge field 5^, which shall be Abelian, not 
spontaneously broken and free from any anomaly. We further assume that the different 
fermions have different charges with respect to this field. The cancellation of the anomaly for 
this gauge field requires that the sum of all the fermionic charges vanishes and therefore 
an anomalous transition leading to the creation of one of each existing fermions perfectly 
respects the 5-gauge symmetry. However, in the case of mixing between the fermions, the 
creation of only light fermions leads to i?-charge violation. 

Note that the mixing between the fermions is only possible if the background has a 
non-vanishing charge with respect to the S-gauge symmetry. It means that the S-gauge is 
broken in the sphaleron or instanton core. This is possible, and is generally the case, even 
if the i?-gauge symmetry is unbroken in vacuum. 

Two points deserve further investigations. Firstly, if the level crossing picture changes 
qualitatively, it is interesting to see if the transition probability undergoes such changes. 
Secondly, we have to understand how charge conservation is ensured. 

These questions are mainly model independent, therefore we choose to resolve them in 
a simple 1+1 dimensional anomalous Abelian Higgs model with two chiral fermions, which 
contains the above paradox. In this particular model, we will show that the probability 
for the creation of two light fermions is zero, unless it is accompanied by the emission of 
some other particles that compensate for the charge asymmetry. The level crossing picture 
must therefore be reinterpreted. We will also see that the nonperturbative transition is more 
probable if there are heavy fermions in the initial or final state, again in contradiction to 
what the level-crossing picture suggests. 

These paradoxes also arise in the electroweak theory. Indeed, the quarks have vari- 
ous electric charges and in the background of the electroweak sphaleron or instanton the 
SU{2) X [/(l)-gauge symmetry is completely broken by the presence of a charged weak field 
background. The resolution of these questions is of great interest for electroweak baryogen- 
esis^. 

This paper is organized as follows. In Sec. m we design a 1+1 dimensional model 
adapted to our purposes. The level crossing picture of this model is derived in Sec. IIIII In 

Even though with the current constraints on the Higgs mass, producing the observed baryonic asymmetry 
within the minimal Standard Model is impossible, electroweak baryogenesis may still work in some of its 
extensions and in supcrsymmetric theories |^. 
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Sec. IIVI we perform a computation of the non-perturbative transition rate in the instanton 
picture and resolve the paradoxes mentioned before. The imphcations of our results for the 
electroweak baryogenesis and a conclusion is given in Sec. |V1 



II. THE MODEL 

We construct here a simple model which contains the paradoxes mentioned in the in- 
troduction. A good candidate for stud ying the creation of massive fermions is the chiral 
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where 1)^0 = - ieA^, 1/(0) = f (|0|^ - v"^) and = Qf) with j = 1,2 labeling the 
different flavors. We use the following representation for the 7-matrices: 



, 75 = 7o7i- (2) 



This model displays similar nonperturbative properties as the electroweak theory. 

We also need a mixing term between the different flavors, which would, at the semi- 
classical level, prevent the energy levels of the fermions from crossing each other. Interaction 
terms between fermions are not possible in this simple model. Indeed, if we write interaction 
terms between different fermions like fij'^\ j^'^ + h.c. we could redefine the fields to 
diagonalize the matrix fij and the theory would inevitably lead to independent fermions^. 
We therefore have to introduce another scalar field, allowing for other Yukawa couplings. 
We will also introduce a U{1)^ gauge field 5^ to give different charges to the two different 
flavors. The bosonic sector reads: 



^ Moreover the first scalar field in vanishes in the center of the instanton, which would allow the levels 
to cross at the center. 
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We have now to specify the charge of each of the four spinor components with respect 
to U{1)^ and f/(l)^. Let us note a]^\ and the charges with respect to and i?^. The 
following choice turns out to serves our aim: 



"iJ ~ 2 ' ~ 2 ' Pi? — 2 ' P. 



and ^' ^' (4) 

o2 _ _e' ^ ' 

L — 2 • 



The gauge symmetries imply that there are two classically conserved electric currents 

fs = + f^R^Rl'^R- 

These currents are in general anomalous but are conserved with our particular choice of 
charges. 

i i 
i i 

The fermionic current 

= W^rn + ^rY'^r- (5) 
is conserved at the classical level, however, its anomaly does not vanish: 

i i 

= ^^.^K- (6) 

This is indeed what we need; there is no gauge anomaly but the fermion number current is 
anomalous. We can now write down an interaction between fermions and scalar field, and 
the fermionic Lagrangian reads: 

The fermionic spectrum consists of two fermions of different mass = vf\ j = 1,2 
interacting with each other by Yukawa coupling to the scalar field x- The vacuum structure 
of the model given by the Lagrangians ^ and Q is non-trivial ^] . Taking the Aq = gauge 



and putting the theory in a spatial box of length L with periodic boundary conditions, one 
finds that there is an infinity of degenerate vacuum states \n), n G Z with the gauge- Higgs 
configurations given by 

^1 = ^, (P = ve'—. (8) 

The transition between two neighboring vacua leads to the creation of two fermions as 
intended: If the vector field undergoes the variation 

^ -1^. (9) 

which corresponds to the difference between two adjacent vacua, the fermionic anomaly is 

SNp = [ ^e^^'F^J^x = ^ / 2dM^x = ^2SAiL = 2. (10) 

J 2lT ^ ZTT J ZTT 

III. LEVEL CROSSING PICTURE 

We build a path in the bosonic field space that goes adiabatically from one vacuum to 
the neighboring one. To this aim, we find the sphaleron and construct a path that relates 
it with the initial and final vacua. Such configurations are relevant for high temperature 
dynamics jl^ . 

Using the Aq = Bq = gauge, the sphaleron in this model reads 



id 



,d 



7VIX 

—ive L tanh(Ma;) 

_ nix , _ 1 / , , \ 

tae L cosh [Mx) 



X 

At = (11) 

Bf = 0, 

with a = \J j^i^v"^ — 2M^). It can be found using results on solitons with two scalar fields 
of Refs. Note that this solution is only valid for a restricted parameter space 

At;^ > 2M2, (12) 
2M2 + Aa^ - /if 2 = 0. (13) 

An example of a path going from vacuum = at r = to vacuum n = — latr = l 
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via the sphaleron at r = 1/2 is 



(jf^ — ve "i"^^ [cos(7rT) + isin(7rT) tanh(Mxsin(7rT))] , 
X^^ — — iae"^^ sin(7rT) cosh~''"(Mxsin(7rT)), 



2nT 



(14) 



Bf - 0. 



These configurations represent a set of static background fields interpolating between vacua 
in which the fermions evolve. The equations of motion for the fermions are 



(15) 



with 



^3 

V *4 = *W 



and 



H 



2 ^ 2 

-idi- 










/20 



v 



(16) 



(17) 



the Dirac Hamiltonian. In the limit of slow transition f ~ 0, the Hamiltonian is time- 
independent and the spectrum of the static Dirac equation i?^„ — En^n for each r leads 
to the level-crossing picture. Of course, an analytic solution to this eigenvalue problem is 
not possible for each r. We therefore give the analytic solutions at a few values of r, check 
with perturbation theory that the interaction potential lifts the degeneracy of the levels 
where they cross each other, and then give the complete level crossing picture resulting from 
numerical computation. 
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A. Fermionic spectrum in the vacuum r = 0, 1 



For r = 0, 1 the fermionic spectrum is the one of non- interacting fermions i = 1,2 



111 ] and is labeled by an integer n. 



K = MlF! + kl (18) 



— r = 



L 

Note that the spectrum is different in the states r = and r = 1. The configuration r = 
is not a true vacuum for fermions, the fermionic contribution to vacuum energy being larger 
for r = than for r = 1. This difference however vanishes in the limit of infinite system 
size lUi]. All states are doubly degenerate in energy except for r = in the case = 0. 



B. Sphaleron configuration, r = 1/2 



H 



(19) 



The Dirac Hamiltonian in the background of the sphaleron reads 

/ -i(9i + ^ iFie-^'rtanh(Mx) -iFg cosh(Mx)"^e"^T^ \ 
-iFie^'r tanh(Mx) idi+^ 

-^^1+^ iFse-^TT tanh(Mx) 

\iF^cos\i{Mx)-^e'^ -zFse^'r tanh(Mx) ^^i + ^ / 

In the F3 = a/3 = case, these equations decouple can be solved separately for each fermion. 
In the limit L ^ 00, one finds two zero-modes, one for each fermion: 



^^0 



e 2L' 



'[cosh(Mx)]^ 
[cosh(Mx)]- 



M 



3 
M 



J = 1,2. 



(20) 



The interaction can be introduced perturbatively. To this aim, the Dirac Hamiltonian is 
separated in two parts H = Hq + W with Hq = H{f3 = 0). In the "^l, subspace, the 
interaction matrix reads 



1 



niTij 



'¥\W\^^) 



il 



with 



nin2 



1 



/ (F3)[cosh(Mx)] 



rL/2 



M 



-1 



dx 




(21) 



and rii = The eigenstates of the matrix Mij are 

= — + \l/2 with energy = I, 
= — with energy = — /. 

We see here that the interaction between the fermions hfts the degeneracy between the states 
and avoids that the levels cross each other. 

C. Numerical results 

The energy levels may be found numerically for each value of r solving the static Dirac 
equation with the Hamiltonian (fT7j) and periodic boundary conditions in the interval of 
length L. 

The results (Fig. |21 show, in the cases of independent and mixed fermions, the creation 
of two fermions (two levels cross the zero-energy line). In the independent case, one of each 
fermion is created (Fig. 12)), whereas two light ones are created in the mixed case (Fig. Oj). 
The latter process violates charge conservation^. For charge conservation to be preserved, the 
transition probability of such a process must vanish. As a precise calculation of the transition 
probability is difficult in the sphaleron picture, we will use the instanton approach in the 
following, which leads to a well-defined semi-classical expansion. Note that the instanton 
picture will be similar to the adiabatic sphaleron transition if the fermionic masses are large 
and their associated time-scale small in comparison to the instanton size. 

IV. INSTANTON PICTURE 

We first derive the Euclidean properties of the model and then compute the transition 
probability for a few representative processes. In Euclidean space, the bosonic Lagrangian 

^ Two light fermions of charge —1/2 with respect to the B gauge field are created 
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FIG. 2: Level crossing of two fermions without mixing, Fi = 0.05, F2 = 0.35, /s = 0, L = 50 and 
h = m = e = e' = l, M = 0.5. One of eacli fermion is created wlien going from r = to r = 1. 



0.4 




0.2 0.4 . 6 . 8 1 



FIG. 3: Level crossing of two fermions with mixing, Fi = 0.05, F2 = 0.35, /a = 0.24, L = 50. Two 
light fermions are created here. 

reads 

+ J 1^1* - ^ i^r + J 1x1' + ^ ixr + ^ ixr - «^). (22) 
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and the fermionic part 



A. Bosonic sector 

In order to find the instanton solution, let us point out the following: if y = i? = 0, we 
know the solution of the remaining equations, it is the Nielsen-Olesen vortex 'l4]. We search 
here for a solution of the same type, adding some generic form for B and x'- 

Xci{r,e) = g{r), (24) 
Bii{r,e) = e'^r>B{r), 

with polar coordinates {it = r = rcos6,x = rsinO), r the unit vector in the direction of 
r and the completely antisymmetric tensor with e^^ = 1. Some details can be found in 
Appendix^ only some results will be given here. An example of profile is given in Fig. 0] 
and the asymptotic form of the different functions are 

/(r) ^ for + 0{r^), A{r) ^ aor + 0{r'), 

9{r) '^90 + 0{r^), B{r) ^ h,r + 0{r'), (25) 

/(r) ''^ 1 + f^K^{V2\vr), A{r) - + a^K,{evr), 

er 

g{r)'^ g^K,{Mrl i?(r) ™ ^, (26) 

where /o, ao, &o, 5'o, /oo, 5'oo, Qoo, &oo are constants found by computing the exact instanton 
profile. 
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FIG. 4: Instanton shape (with a different scale for the field B) for the following values for dimension- 
less constants (see Appendix^l): m'^ = ^ = 1, /i=^ = 4, ^' = ^=4, p = ^ = 1, H = j = 3. 



B. Fermions 



The fermionic fluctuations in the background of the instanton (jz^ are given by if \1/ 
with: 



H 



Hi 



Ho 



Hi h 
1 1 H^^ 

-tfir tdo -di + l i-Ao - tAi) + |(5o + iBi] 

-tdo -di + l i-Ao + lAi) + I {-Bo + tBi) iM 

-iM* tdo -di + +l i-Ao - tAi) + I {-Bo - iBi) 

-zdo -di + l i-Ao + tAi) + |(5o - tBi) 



A^l 1.4' " I. (27, 

y \ z/sx 

The zero-modes are found solving the equations = 0, with H the Dirac operator in 
the background of the instanton. We use polar coordinates (r, 6) and expand fermionic 
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fluctuations in partial waves \& = ^^=-00^*^^™- "^^^^ leads to the following equations: 

0, 



dr r 



d m — 1 
dr r 



0, 
0, 



In the case /s = and B{r) = 0, the two fermions decouple and their zero modes are 



ip^{r) oc 



exp 



- / dr'{vf,f{r) + -A{r)) 



, J = 1,2. 



(28) 



(29) 



If /s 7^ the zero modes have the following asymptotics^ 



-Fir 



\ 



Fzr^-ie 



-Pie 



012 



J 



\ 



(30) 



/ 



where ai and 02 are normalization constants and f3i and vanishes in the decoupled /s = 
case (j2ni)- In the coupled system, the existence of the zero- modes can be checked using the 
method of Ref. but /3i and (32 have to be computed numerically. The coefficients /5i,2, 
which will be needed to compare transition probabilities, are given as a function of the 
fermion masses in Fig. El Within some approximations, the constants /3i,2 can be found 
analytically (see Appendix . 



C. Transition probability 

We start with two decoupled fermions (/s = 0) and introduce the interaction pertur- 
batively. It is clear what happens here; the interaction term —if^'^^^-^'^'^x + h.c. allows 
for the decay of the heavy fermion into a light fermion and a x boson, a process which 
conserves the charge and which can be taken into account in final state corrections (see 

^ We consider here the approximation B = (or e' = 0), which does not lead to observable changes (see 
Fig. EI) 
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FIG. 5: Coefficients (32 (dashed lines) and Pi (triangles) as a function of the mass F2 for some 
different light fermion masses -Fi = 0.01, 0.1, 0.5, 1 and for = 0.2. The line represent Pi = P2 
in the degenerate case Fi = F2. The constants -Fi,2) /s are in units of ^/Xv. 

[ifij ]. for inclusion of fermions see [itI). This is not what we we are interested in here. In 
the non-perturbative regime, two light fermions are created in the x background, where the 
U{1)^ gauge is broken and a x boson should be emitted from the instanton tail as the U{1)^ 
gauge symmetry is restored far from the instanton center. We will show here that processes 
violating charge conservation have vanishing probability. 



1. Green's function 

Green's functions with creation of two fermions and an arbitrary number of other particles 
read 

G''\xi,X2, yi-.-Vn, Zl...Zm, ^1, W«) = (31) 

/n m I 

v^v^vxvx*vr^ e-^[*'^'>^'>^*'^] ^''{xl)^\x2)l[x{y^)l[^*{yJ)l[v{yk), 
i=l j=l k=l 

where rj stands for all neutral bosonic degrees of freedom, nfc"=i viy'D contain the field 
A, (f) and neutral pairs of fermions and the variable is a spinor containing the two fermions 
as in (fTH|) and a,b = 1, .A. 
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FIG. 6: Coefficients /3i^2 as a function of tlie mass F2 for some different couplings /s and for 
/i = 0.01. Tlie constants -Fi,2) /s ai^e in units of \/\v. 

The main contribution to the Green's function for the creation of two fermions comes 
from the sector with one instanton (g = —1)^. In this sector, fermions have two zero-modes 
flH()|l . The Gaussian path integral over fermionic degrees of freedom can be evaluated, leading 
to the fermionic determinant with zero-modes excluded and the product of the fermionic 
zero-mode wave functions^, 

G{xi,X2,yi...yn,zi...Zm,wi,...,wi)= VxVx*Vr] (32) 

Jq=~l 

n m I 

i=l j=l k=l 

2. Collective coordinates in the one instanton sector 

The bosonic action is expanded around the instanton configuration. Gaussian integra- 
tion over the quadratic fluctuations gives a determinant det^Df^J^^ . However, zero modes 
associated to symmetries require introduction of collective coordinates. There are two trans- 
lation zero modes and one coming from U{1)^ gauge. Performing an infinitesimal global 

^ More precisely, in the dilute instanton gas approximation the result from the one instanton sector can be 

exponentiated to give the complete Green's function [l^ . 
® Note that the zero- mode functions still depend on the background = i^'^lx, X*tV]- 
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gauge transformation, we get 

5x = e'^X-X-^Pe''g{r). 5<P = 5A = 5B = Q. (33) 

Note that the Uil)"^ gauge is broken, there is no normahzable zero mode associated to this 
symmetry. Rotation symmetry does not lead to a further zero mode.^ 

Collective coordinates are introduced as follows. The integral over the translation zero 
modes are replaced by an integral over the instanton position xq. The integral over the U{1)^ 
gauge zero-mode is replaced by an integral over all possible global gauge transformations (3. 
The Green's function reads: 

»,.....«,,) = //.o dP det'(A-„„,)iV«iV,. (34) 

n m I 

xdet'(L'D"5^J(a;i - a;o)^/'J(x2 - xq) JJ x^iVi - xq) JJ x*(i{yj - xo) JJ v^iVk - xq), 

i=l j=l k=l 

with 

X,3 = e*^Xd, Xl = e-'f'xci, (35) 
Vf3 = Vci, V'^ = e^^^^^^„ j = l,2, (36) 

and Nb, Nfr the normalization factor coming from variable change to collective coordinates. 
To simplify the notations, we also introduced the matrices Fj, i = 1, 2, 5 acting on the four 
dimensional spinor (fT^ as: 

r.J''°], r.J""), r.^f"'^ ° ], (3T) 

\^ y \o hj \o -hj 

where I2 is the identity on a two dimensional subspace. 
3. Fourier transformation of the Green's function 

The Fourier transformation of the Green's function after integration over the instanton 
location xq reads (writing spinor indices explicitly) 

G''\h,k2,pi...pn,p[...Pm,qi,-,qi) = (2vr)V') (P) / dPK (38) 



X 



i=l j=l k=l 



Rotations give the same zero- mode as U{1)^ gauge transformations 
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where k = e"^- det'(if,„,i)iVBiVi,det(DL)"^ and P = h + k2 + Y.l=lPi + YZlP'^ + Y!^=l<l^■ 
The integration over the instanton location leads to momentum conservation. In a similar 
way, integration over gauge rotation (3 enforces charge conservation. Indeed, the integral 
over (3 is non-zero only if the powers of e^^ cancel, that is to say, if charge with respect to 
the gauge field is conserved. 

As the different components of the spinors have different powers of e^^ , different cases 
have to be considered. We will concentrate here on three interesting situations, from which 
we will be able to derive some general conclusions. 

D. Examples of allowed matrix elements 

First consider a process involving one scalar as initial state, which decays into two 
fermions. In this case the integration over the coordinate (3 leads to: 

G''\h, h, qi) = (27r)25(2) (P) «:0,,(gi) (39) 
X {(v,ijl,{k,)Y (T,ijl{k,))\ (T,^l,{k,)Y (V^i^lik,))'^ . 

Applying the reduction formula, we get a non-vanishing matrix element for two different 
fermions only by multiplying the Green's function by two fermionic legs v}{ki), u^{ki), 

tM{kl, k2, ql) = (27r)25(2) {q - k^ - k^) k t{q^ + m]j)^a{q) (40) 
( iu\k^){k^ + F)T^^l,{k^) tu\k,){h + F)r,^l{k2)) . 

A straightforward calculation gives (see Appendix E)) 

|M(M,A;2,gl)|2= (2(27r)Vocaia2(l+/3i/32))'. (41) 
The decay rate is after integration of the phase space (supposing nii << m^): 

^'t' = V- f dUps\M{kl,k2,ql)\^ = - -4 ^(2(27r)Voo«i«2(l + /?i/?2))'. (42) 

Secondly, we consider a process involving one x scalar as initial state. The Fourier 
transformation of the Green's function reads 

G'^'ik,, q) = (27r)25(2) {q-k,- k,) k xlM) [^li^liihif (Fi^^ (fc^))' ■ 

-"^^ Note that this do not depend on the existence of the field, but on the existence of the associated global 
symmetry. Therefore the requirement of charge conservation will persist in the limit e' — > 0. 
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Applying the reduction formula, we get the matrix element for the creation of two light 
fermion by multiplying the Green's function by two light fermion legs u^{ki), u^{k2), 

tM{kl, k2, ql) = (27r)25(2) (g - fci - k^) k t{q^ + ml)x*ci{q) (43) 
iu\ki){h + F)T,^li{k,) iu\k2){k2 + F)T^i;l{k2). 

A straightforward calculation gives (see Appendix [nj) 

|M(fcl, fc2, gl) p = {2{2^fKg^aia2(32f. (44) 

The decay rate after integration of the phase space (supposing nii « m^) reads 

{2{2n)^ Kg^aia2l32f 
Tx = / ^ ■ (45) 



A similar process involves the scalar x*, which decays into two heavy fermions: 

^ ^ {2{27rfKgMa2(3if ^^^^ 
2m'^^m'^ — Am\ 

Generalizing these three examples, we see that any interaction leading to the creation of 
two light respectively heavy fermions contains a factor (32 respectively (3i. Processes leading 
to the creation of one of each fermion contain a factor 1 + /?i/?2 ~ 1- 



E. Discussion of the different transition probabilities 

The integration over the collective coordinate associated with the gauge symmetry U{1)^ 
leads necessarily to charge conservation. Therefore the process described by the level cross- 
ing picture (Fig. ^ cannot take place without the emission of some other particle that 
compensates the additional U{1)^ charge. The possible initial and final states are more 
restricted than suggested by the level crossing picture. 

We shall now discuss the transition probability of allowed processes. We leave aside for the 
moment the phase space factors, they are not large in 1+1 dimensions. The main factors 
that distinguish the transition rates (j42p . (j45j) . (j46j) for the three possible fermionic final 
states are the constants /3i^2- This is also true if more complicated processes are considered. 
As expected, if the fermions are light and weakly coupled, the probability to create one of 
each fermion is much larger; it is proportional to 1 + /3i/52 = 1, see Fig. El However, in 
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the case where one fermion is very heavy, it is unexpectedly favored to create two heavy 
fermions. The creation of two hght fermions suggested by the level crossing picture is indeed 
suppressed. In the case of slow transitions (heavy fermion masses Fi 2 or large instanton 
radius Vinst) the probability of creating light fermions is larger, but it is still suppressed if 
the mass hierarchy is large (-^ ^ 1). 

V. CONCLUSION 

In the model considered here, the level crossing picture suggests a particular transition 
which must not and does not occur. A possible way out would be to reinterpret it as 
follows. The level crossing picture only knows about fermions and the correct bosonic 
content of the initial and final states should be added by hand when dealing with a physical 
transition. More precisely, all symmetries that are broken by the fermionic initial and final 
states should be restored by supplementary bosonic operators. However, even with this extra 
requirement, the level crossing picture suggests the creation of two light fermions, a transition 
that turns out to be suppressed. Furthermore the most probable transition, computed with 
the inclusion of the first quantum corrections, would impose the energy levels to cross each 
other several times on the level crossing picture in spite of the interaction potential. Note 
that this is perfectly possible in quantum field theory although forbidden in the adiabatic 
quantum mechanical description. 

The results for the transition probability are rather surprising; for heavy fermions, such 
as the top quark, or adiabatic process rinstF2 ^ 1 (Sphaleron at high temperature) the 
probability of creating two heavy fermions is large. In the realistic electroweak theory, the 
phase space factor may be dominant and may change this conclusion. It is therefore very 
interesting to reproduce similar computations in the frame of the electroweak theory at high 
temperature, or at high energies. 

A more interesting setup would be to include heavy quarks in the initial states. The phase 
space factor as well as the matrix element are then large. In this case, the nonperturbative 
transition rate can be enhanced by a huge factor (see Fig. Ej). A high top quark density 
could therefore catalyze the nonperturbative transition rate. This phenomenon is relevant for 
baryogenesis at the electroweak phase transition. It could provide a mechanism to enhance 
the baryon number violating transition rate in the symmetric phase, while suppressing it 
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in the broken phase. Indeed, while bubbles of true asymmetric vacuum expand in the 
symmetric universe, it may be that top quarks are more reflected by the bubble wall and 
are rare inside the bubble, and over-dense outside. This density asymmetry will render the 
nonperturbative rate faster outside the bubble, while slower inside. 

It should be noted that the present calculation deals with the instanton rate, although 
at high temperature, the sphaleron rate is the relevant quantity. It would therefore be very 
interesting to flnd out if the sphaleron rate also displays these interesting features. 
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APPENDIX A: THE INSTANTON 

The bosonic Lagrangian in the d^^A^ = d^B^ = gauge gives the following equations 
of motion, 

- d^d^4> + 2ieA^d^4> + e^Al<p - Xv^4> + X \4>f <j) + h = 0, 

-d^d^X + 2^(eA^ + e'B^)d^x + {eA^ + e'B^fx, 

+M\ + K\x\^X + hm^ -V^)X = 0, 
-d^d^A, + + X* V^X) + e'A^(|0|' + Ixl') + ee'B^ = 0, (Al) 

-d^d^B^ + i''-X*^^X + e'^B^\x\^ + e'eA^\x\^ = 0. 

We are looking for a solution of the type ()24|) . As for the Nielsen-Olesen vortex, we impose 
the asymptotic behavior of the functions A and /: 

/(r)^/ir, f{ry^l, A{r) a,r A{ry^ ^. (A2) 

er 

For the flniteness of the action, the function g(r), B(r) should respect the following boundary 
conditions: 

i?(r)^6ir, g{r)'^g,, B{r) + rB\r)'^ ^7(r) ™ 0. (A3) 
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We also introduce dimensionless variables with the substitutions 



A = A 



7VX 



f = f 



9 



e e 
The remaining parameters are 

A , A 



5, B = B— 



A h ^ 



The equations of motion ()A2|) in polar coordinates and with the ansatz 

2 



- fir) - I fir) + 



A{r) 



fir) 



+/i/>)'-/>)+^W7(r) = 0, 



1 



-9"ir) - -g'ir) + {A{r) + B{r))'g{r) + pg{r)' 
r 



1 



A{r) - - 



+ —[B{r) + A{r)]~g\r) = 0, 



reads 



(A4) 



(A5) 



(A6) 



where the prime means derivative with respect to r. 



APPENDIX B: ANALYTICAL APPROXIMATIONS FOR THE FERMIONIC 
ZERO-MODES 

For small coupling /a, and small instanton size a in units of fermion mass, we can get a 
rough approximation by perturbation theory. We checked numerically, that it corresponds 
reasonably well to the exact case and will be sufficient for the following discussion. We are 
interested in the case were the first fermion is very light in comparison to the second one 
and in comparison to the scalar field, Fi <^ m^. 

The first step is to eliminate the field by the variable change \1/ — > exp (— | / drA{r)) \l/ 
and contract the four first order differential equations into two second order ones. One 
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obtains the following equations for the new variable 



m(m — 1) f'(r)m „, ,0 /.2\ t / x /'('")^2('^) t/// x 



r/(r) y /(r) 

= fir)gir)fj,^,ir), (Bl) 
m(m-l) f(r)(m-l) ,,,2A.T.r^ /'W^sW 



r/(r) 



/(r)V2^) Mr) - ^'^"j^f^ + mr) 

f{r)g{r)hh^,{r), (B2) 



1 f^,, m\l/2(r)\ ^ , , 1 /(m - l)\l>3(rl 



with 

We discuss the case of the zero- mode in m = partial wave first. At zero-order of 
perturbation we have the two first components {'^'12) given by ()29p and the two last ones 
(V'34) vanish. If we consider now a non- vanishing f^g{r) in ()B2|) . the function ip^ is given at 
first order perturbation theory by 



^3(r) = /2/3 j dr'G{r,r')f{r')g{r')ijl{r'), (B4) 

where G{r,r') is the Green's function of the differential operator in the left hand side of 
equation ()B2jl . We were not able to find a general expression for G{r,r') for an arbitrary 
function /(r), but satisfactory results are obtained using the Green's function G{r,r') for 
constant^^ /(r) = v. In this particular case 

G{r, r') = sinh(/2r<) exp(-/2r>), (B5) 

J2 

with r> = maxfr, r'), = min(r, r'). Form (IB4I). we get, for r ^ 1: 



^3(r) = -/3t^exp(-/2r) J dr' smh{f2r')g{r')^,{r'), (B6) 

If the inverse fermion mass is small in comparison to the typical extent Vinst of the 
function g{r), we have: 

roc roc 

Pi = f3V dxg{x) smh{f2x)e~^^'' fsv dxg{x) sinh{f2x) (B7) 
^0 Jo 



This approximation is exact in the hmit of small instanton size and precise for light fermions, because 
they do not probe the instanton center. 
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In the case of a large mass /2 and large instanton size, the constant jSi can be large from 
the presence of the sinh. In the case of small /2, the integral can be further simplified to 
A = fzF2 J xg{x)dx. A similar computation^^ can be performed for '^p2^ 

= f^v j dx g{x) sinh(/ia;)e--^2^. (B8) 

If fifinst ^ 1 we have /32 = /3-F1 / dxxg{x)e~f^^ , which may be large for a large instanton 
size, if it is not suppressed by a large fermion mass /2. 



APPENDIX C: FOURIER TRANSFORMS 



For computing cross sections, the unitary gauge is best suited. It is however known to 
be singular, which may lead to discontinuities in the fermionic wave functions. This can be 
easily cured using the following regularized gauge condition: 



a{r,9) =e-27re'{9-7r), 



(CI) 



where Q'^{9 — n) is continuous and goes to the step function as e (see Ref. fore 
more details). The Fourier transforms^^ of the fields are (we consider here only the case 
/s = for fermions): 

foo 



(f){p) = 
Mp) 



,2 ' 



xip) 



9c 



m 



p^ + 



(C2) 
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F,+p 



TCp 



Ff+p^ 



where p = ^/p^p^ and 9g the angle between the spacial axis and the vector p. If /a 7^ but 
the field B is neglected, the Fourier transforms of the two fermionic zero modes read: 



i'liip) 



—ia\ 
—iai 
—ia2 
—ia2 



e„ Fi+p \ 



e 2 



9r, F-z+p 

2 F2+P 

F^ ) 



( 



^Up) 



\ 



-t(3i 



np F?+p^ 
2 



\ 



F^ 

^^^^e,F^ 

Trp 



(C3) 



F^+p2 y 



The relevant green's function is given by l|B5(l . were we replaced /i — ^ /2- 
13 We retain only the pole term here, the rest do not contribute to the final amplitude [l9j 
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